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Abstract. 

We study experimentally and theoretically the sedimentation of gels made of 
strongly aggregated colloidal particles, focussing on the long time behavior, when 
mechanical equilibrium is asymptotically reached. The asymptotic gel height is 
found to vary linearly with the initial height, a finding in stark contrast with a 
recent study on similar gels [Manley et al. 2005 Phys. Rev. Lett. 94 218302]. We 
show that the asymptotic compaction results from the balance between gravity 
pull, network elasticity, and solid friction between the gel and the container 
walls. Based on these ingredients, we propose a simple model to account for the 
dependence of the height loss on the initial height and volume fraction. As a result 
of our analysis, we show that the static friction coefficient between the gel and 
the container walls strongly depends on volume fraction: the higher the volume 
fraction, the weaker the solid friction. This nonintuitive behavior is explained 
using simple scaling arguments. 
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1. Introduction 

Most colloidal systems tend to aggregate, unless special care is taken to stabilize them, 
e.g. by designing particles with surface charges or a sterically stabilizing polymer layer 
that prevent them to approach too closely, where short-ranged, strong Van der Waals 
or hydrophobic/hydrophilic attractive interactions prevail. As a result of aggregation, 
ramified structures, often with a fractal morphology, are formed. For high enough 
particle volume fraction, tp , and interaction strength, these aggregates form a space- 
filling network which is termed a colloidal gel 1 . The fate of colloidal gels depends 
crucially on gravity, since the magnitude of both their clastic modulus and yield stress 
are generally comparable to the gravitational stress due to the mismatch between the 
density of the particles and that of the solvent in which they are dispersed. 

Because colloidal aggregates and gels are ubiquitous in the food, drug, cosmetic 
and oil industry and in biological systems, a large number of works has been 
devoted to the sedimentation (or the creaming) of aggregated colloidal systems. 

§ To whom correspondence should be addressed (lucacip@lcvn.univ-montp2.fr) 
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Most investigations have dealt with "weak" gels, where the interparticle attraction 
is comparable to k^T, the thermal energy El El- Experimentally, this is typically 
the situation encountered in systems where the attraction between the colloids is 
induced by the addition of smaller particles or polymers via the so-called depletion 
mechanism @. For these systems, various sedimentation regimes have been observed, 
ranging from "creeping" sedimentation (a linear decrease of the gel height, h, with time 
t) to "delayed" sedimentation (h(t) is nearly constant during an initial latency time, 
followed by a sudden collapse of the gel and a final slow compaction of the sediment) . 
Moreover, the sedimentation of samples with identical composition may drastically 
differ according to the container size and shape and in particulary depending on the 
initial height of the gel, hg. This behavior can be rationalized according to whether 
or not stress can be transmitted through the sample over distances comparable to 
the container size [3 El- This is the case for the strongest gels, where the attractive 
interactions are large enough to make it likely that connecting paths along the gel 
network persist, in spite of the continuous breaking and reforming of particle bonds 
due to thermal energy. 

The group of C. Allain has extensively studied stronger gels, where particle bonds 
are unlikely to break because of thermal energy |§1 1101 ITT1 112) . These gels, however, 
are still weak enough to experience extensive restructuring due to the gravitational 
stress. For the most diluted samples, the very formation of a gel is prevented, because 
clusters grow more slowly than in concentrated samples and they sediment before 
getting a chance to form a system-spanning network. By contrast, Manley and 
coworkers have studied "strong" gels made of silica particles where neither thermal 
energy nor gravitational stress are large enough to cause significant restructuring of the 
network They find that h(t) relaxes exponentially towards an asymptotic height 
hoc and propose that this behavior results from the competition between network 
elasticity, gravitational pull, and viscous drag of the solvent through the gel pores. 

In this paper, we study strongly attractive colloidal gels similar to those of Manley 
et ai, albeit at slightly smaller volume fractions. Although we observe a similar 
exponential evolution of h(t), we find that Aoc = Iiq — hoc is proportional to the initial 
height, a behavior incompatible with that predicted by the model of reference [T3|. The 
volume fraction dependence of the sedimentation is also found to deviate with respect 
to that postulated in 13 . We propose a more general model for the equilibrium 
behavior of strongly attractive colloidal gels: in addition to network elasticity and 
gravitational stress, our model introduces a term describing the solid friction between 
the gel and the cell walls, using a formalism similar to that first developed by Janssen 
for granular materials The model accounts for all our experimental observations 
and allows us to estimate the static friction coefficient, fj, s , between the gel and the 
walls. Quite surprisingly, we find that fx a decreases with increasing volume fraction. 
We propose a simple scaling argument to show that this counterintuitive behavior is 
due to the fractal nature of the aggregates forming the gel. 

2. Materials and methods 

We use aqueous dispersions of charge-stabilized silica spheres (Ludox TM-50) with 
radius a « 11 nm. Aggregation is induced by mixing equal volumes of a particle 
suspension and a NaCl solution to a final salt concentration of 2M and final particle 
volume fractions ipo ranging from 2.5 x 10 -3 to 10~ 2 . At this ionic strength, 
electrostatic repulsion becomes negligible and the particles experience short-ranged, 
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strong van der Waals interactions, leading to the formation of aggregates in the 
diffusion limited cluster aggregation (DLCA) regime JS]- The bond energy is of the 
order of several tens of k^T, making thermally activated rearrangements extremely 
unlikely. The aggregates have a fractal dimension df w 1.9; during aggregation 
their typical size increases until they form a space-filling network. When this gelled 
structure is formed, the aggregates have a size distribution peaked about an average 
cluster size £ ps aipQ^ df ^ The gelation time can be calculated within the DLCA 
model: for our gels it is less than a few seconds |17j . 

The gels are viscoelastic materials: oscillatory strain measurements show that 
over a large range of frequencies (typically 10~ 2 rad/sec < ui < 10 rad/sec) the 
storage modulus G' is nearly frequent-independent and about one order of magnitude 
larger than the loss modulus G" |201 12 1| . Due to the fractal morphology of the gels, 
G' increases very strongly with volume fraction: G' oc p^ where v 3.6 — 4.0 
|18l 1191 I2UI I21j . Moreover, the gels are able to sustain large shear strains up to 
about 10% before the network is significantly altered and linear elasticity fails |20j. 
Applying larger strains leads initially to strain hardening, but eventually results in the 
disruption of the gels, although gel breaking may not occur instantaneously |20l I13j . 

The particle density is 2.35 g/cm 3 , while that of the solvent is 1.15 g/cm 3 (due 
to the high salt concentration). The density mismatch, Ap = 1.2 g/cm 3 , is large 
enough to make the gels sensitive to gravitational stress, leading to a slow compaction. 
Experiments are performed by directly preparing the gels in glass cylindrical cells of 
diameter D — 1.7 cm and letting them settle unperturbed. Note that contrary to the 
protocol of reference |10j no stirring is applied after the initial mixing of the particle 
and salt solutions. The gels are imaged with a CCD camera to capture the settling 
kinetics. We check that no cracks appear in the bulk of the gels. As observed in other 
experiments [30, at the earliest stages of sedimentation the meniscus at the top of the 
suspension is emptied from the particles. Once only solvent is left in the meniscus, the 
gel starts settling uniformly throughout the cross section of the cell. A reasonably flat, 
sharp interface is observed between the highly turbid gel and the clear supernatant, 
allowing us to measure the time evolution of the height of the gel, h(t), with an 
accuracy of about 5%. We define t = as the time when the meniscus is emptied and 
the uniform settling starts. We study gels at six volume fractions: (po = 2.5 x 10 -3 , 
3 x 10~ 3 , 5 x 10~ 3 , 6 x 10~ 3 , 7.5 x 10" 3 , and 10~ 2 . For p < 2.5 x 10~ 3 the gels often 
develop cracks and the experiments are less reproducible, while for po > 10 -2 the gels 
hardly settle, making a precise determination of h(t) impossible with our apparatus. 
For each ipo, we vary the initial height of the gel, ho, in the range 0.5 cm < ho < 5 cm. 

3. Experimental results 

Figurenshows the evolution of the gel height measured for ipo = 0.3 x 10~ 2 and various 
initial heights ho- The time evolution of h(t) is well fitted by a simple exponential 
relaxation: 

ha - h(t) = Aco [1 - expH/r)] (1) 

where A^ — \hoo — ho\ denotes the asymptotic settling and /i^ = hit — » oo). A 
similar behavior has been reported recently for the "strong" gels of reference |T3"| . 
Both the asymptotic settling and the characteristic time r grow with the initial height 
of the gel; for the tallest gels (h = 4.9 cm), r is as large as a few hours. The same 
exponential relaxation, Eq. Q, describes the height evolution of all the gels that we 




Figure 1. Time dependence of the settling of a gel with ipo = 0.3 X 10 — 2 . From 
top to bottom, the initial height of the gel, ho, is 1.23, 2.45, 3.68, and 4.9 cm. 
The lines are fits of Eq. Q to the data. 
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Figure 2. Asymptotic settling of the gel height Aoc (a), and characteristic time 
t for the sedimentation (b) as a function of the initial height ho- From top to 
bottom, <p = 2.5 x 10" 3 , 3 x 10" 3 , 5 x 10" 3 , 6 x 10" 3 , 7.5 x 10" 3 , and 1 x 10" 2 . 
Both quantities vary linearly with ho, as shown by the fits (straight lines). 



have studied, except for some samples where fractures appeared during the settling. 
For fractured gels, h(t) deviates from an exponential relaxation and the asymptotic 
settling is larger than that expected for pristine samples. In the following we analyze 
only data for non-fractured gels. 

We plot in Fig. |2 the h dependence of the asymptotic settling (a) and of the 
characteristic settling time (b). Remarkably, both Aoo and r vary linearly with the 
initial height of the gel, as shown by the lines that are fits of form = aho + a.\ 
and r = (3h + 0\, respectively. As one can see in Figure |U both ol\ and /3i are almost 




Figure 3. Left axis, solid squares: tpo dependence of the proportionality 
coefficient a between Aoo and ho extracted from the linear fits of Fig- Et) • Th e 
line is a power law fit yielding an exponent — f. 54 ±0.17. Right axis, open circles: 
proportionality coefficient j3 between r and ho extracted from the linear fits of 
Fig. Eb). The line is a power law fit yielding an exponent —1.62 ± 0.11. Inset: 
ratio a//3 vs ipo. Note the weak variation of a//3 with volume fraction, indicating 
that the initial sedimentation velocity is almost constant, as discussed in the text. 



zero, making the fits indistinguishable from straight lines through the origin. In fact, 
we find that (3\ = is compatible with experimental uncertainty, while ai is slightly 
but consistently smaller than 0, a feature that will be commented further below. 

For a fixed cell diameter and material, the prefactors a and (3 depend solely on 
volume fraction. Their tpo dependence is shown in Fig. [21 in a log- log scale. The data 
can be fitted reasonably well by power laws: a ~ tpy and (3 ~ <pg, with u = — 1.54±0.17 
and v = —1.62 ± 0.11. Quite intriguingly, u m v; this indicates that for t — > the 
settling velocity of the gels is nearly independent of both their height and volume 
fraction. Indeed, for t << r Eq. Q yields \dh/dt\ w Aqo/t = a//3, which varies 
weakly between 1.45 x 10~ 4 and 2 x 10~ 4 cm/sec as shown in the inset of Fig. [3] 

In order to understand the behavior of our gels, we start by comparing the results 
presented in Figures ^ El with previous works. C. Allain and coworkers have measured 
the quantity fi = 1 — A^/Zip = 1 — a and found that in the so-called "friction- 
dominated" regime f2 is independent of the initial height ho [§1E11> m agreement with 
our observation Aoo ~ ahg. However, they find Q oc ipo, a scaling that is incompatible 
with the power law a ~ (p~ 1S2 shown in Figure|31 as we also check by plotting directly 
fl vs ipo (data not shown). This discrepancy suggests that the physical mechanism 
leading to the settling of our gels is different from that of references [HIEI- Indeed, the 
gels studied by Allain's group experience extensive restructuring and fractures during 
sedimentation |1UI12| . in contrast to our gels. During sedimentation, the gel does not 
retain its integrity and the scaling £1 oc ipo results from the formation of a sediment 
composed of randomly packed aggregates issued from the gravity-induced breaking of 
the gel. 
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An exponential relaxation similar to that shown in Figure ^ has been reported 
for gels made of silica particles by Manley et. al, who studied the regime where no 
fractures are observed |13| . similarly to our experiments. This sedimentation kinetics is 
interpreted as the result of the interplay between gravitational pull, network elasticity, 
and viscous drag of the solvent through the network pores. However, our results differ 
distinctly from those of reference ^Sj in that we find i) Aoo oc ho rather than Aoo oc /i 2 , 
as predicted by the model of ^Hl, h) t oc ho, as opposed to r oc h"o, in) r oc y^ 1 ' 62 , 
rather than r oc (1 — (po); iv) dh/dt\t=o ~ constant, while Manley and coworkers find 
dh/dt\ t= o oc <^( 1 ~ <i f)/( 3 ~ d f) rj ip~ - 8 . Therefore, it is clear that in the case of our gels 
a fundamental physical ingredient is missing in the model of reference |13| . 

4. Model and comparison with the experiments 

We propose that, in addition to the mechanisms listed above, solid friction between the 
gel and the cell walls also plays a crucial role in the settling of our samples. A hint of 
the importance of solid friction is provided by the observation that both and r are 
modified by using plexiglas cells rather than glass cells [22J. Moreover, the importance 
of solid friction is consistent with the arguments of reference [H] , where it was argued 
that, for gels with strong interparticle interactions, the characteristic length over which 
stress can be propagated exceeds the container dimensions. Accordingly, "strong" gels 
such as those studied here should feel the influence of the cell walls. Solid friction is 
particularly appealing since it provides a simple physical mechanism to explain the 
linear scaling of A m with initial height ho, as shown by the following qualitative 
argument. Solid friction screens the gravitational load acting on any given cross 
sectional plane S in the gel: because part of the weight of the gel column laying above 
£ is sustained by the walls, it is as if in E the gel experienced only the weight of a 
column segment of effective height L, shorter than the whole column. One can then 
imagine to divide a gel of height ho into segments of height L: under gravitational 
stress, each segment is compressed by the same amount, regardless of the number of 
segments that compose the gel. Therefore, the total compression is simply proportional 
to the number of segments contained in the gel, i.e. Aoo oc ho/L. 

To make this argument more quantitative, we build a model for the settling of the 
gels, assuming that they are elastic media that compresses under their own weight, 
without any restructuring, and subject to solid friction against the cell walls. The 
rate of compression is limited both by the viscous friction of the fluid through the 
network, which can be viewed as a porous medium I13| . and by solid friction. In 
the following, we will focus on the behavior for t >> r, when mechanical equilibrium 
is asymptotically reached and gravitational stress is counterbalanced by both the gel 
elasticity and solid friction; the kinetics of sedimentation will be addressed in future 
work. We choose the z axis along the vertical direction, with z = the bottom of the 
gel and z — ho its top at time t = 0. Due to gravity, the gel is submitted to a uniaxial 
compression p(z) = — a zz where o~ zz is the vertical stress. 

We first recall the behavior in the absence of wall friction. Newton's law for a gel 
slice of thickness dz at height z yields 

-Apg<p(z)—dz- 1 f—dz = 0, 2 
4 az 4 

where tp(z) is the local volume fraction, D the cell diameter and g the acceleration of 
gravity. The first term describes the gravity pull exerted on the slice, while the second 
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term accounts for the elastic response of the material. The conservation of the total 
number of particles imposes the additional condition 

rhoo rho 

I (p(z)dz = / cp dz =h Q ip , (3) 
Jo Jo 

while the boundary condition is 

<p(h oo )=<p , (4) 

since the top layer of the gel is uncompressed. If the deformation is not too large, the 
pressure and the volume fraction are related by the equation of state of an elastic 
medium in the linear regime: p(z) — C tp ^~ ipo , where the uniaxial compression 
modulus C is a linear combination of B, the bulk modulus, and the shear modulus: 
C = B + AG/3. Equation J3J) can be solved for ip(z); by imposing the boundary 
condition, Eq. l|4"|l. one obtains the asymptotic profile: 

<p(z) = <^ exp [(/ioo - z)/X] , (5) 

where we have introduced A = C/(Apg(po), the characteristic length over which an 
elastic material is deformed due to gravity (the stiffer the gel, or the weaker the density 
mismatch, the larger A). By requiring that Eq. (JSJ) satisfies particle conservation 
(Eq. |3J), one finds Aoo = ho — Aln(l + ho/X). In the limit ho << A, this relation 
simplifies to Aoo = hf ) /(2X) + XO(ho/X) 3 , whose leading term coincides with the 
expression given in reference We estimate G ~ 7 Pa for our gel at cpo — 1 x 10~ 2 
by extrapolating the data of ref. [21] • By assuming that the order of magnitude of C 
is the same as that of the shear modulus plateau, one estimates A « 13 cm. [Using 
the full model that takes into account solid friction, we will obtain A directly from the 
experimental data, finding a result consistent with this estimate, see Fig. [S] below]. 
Hence, ho < X in all our experiments (indeed, in most cases ho << X) and a quadratic 
scaling of A m with ho should be observed. By contrast, as recalled above, we find 
Aqo oc ho, indicating that gravity and elasticity alone are not sufficient to correctly 
describe the behavior of our samples. 

Solid friction has been first modeled for a container filled with a granular material 
by Janssen in 1895 14 . In the Janssen model, the pressure p(z) is partially redirected 
in the horizontal direction leading to an horizontal stress o~ rr (z) = Kp{z) with 
K = g + 4g/3 the redirection coefficient. When mechanical equilibrium is reached, 
the effective static friction of a slice of material of thickness dz confined in a cylinder 
of diameter D is oriented upward and has magnitude irfi s DKp(z)dz, where /i s is the 
static friction coefficient between the grains and the container wall. The Janssen model 
predicts that, when moving from the top of the grain column downward, the pressure 
increases reaching exponentially a saturation value with a characteristic screening 
length L = D/(Afi s K). This is very different from the case of a liquid where the 
pressure increases linearly with depth: the saturation is due to the fact that the 
cylinder wall supports part of the weight of grains: the contribution to the vertical 
stress on a horizontal plane E due to the material at a height Az above E decreases as 
exp(— Az/L). Hence, the larger the contribution of solid friction, the smaller L. The 
universality of the Janssen scaling for the stress saturation curve has been recently 
shown by Ovarlez et al. [231 - Moreover, the Janssen model has been generalized to 
dynamical situations |24| . 
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We propose to describe the solid friction between the gels and the cell walls using 
Janssen's formalism. At equilibrium, Eq. @ becomes 

ttD 2 dpnD 2 

— Apgip(z) dz — — dz + ir[i s DKp(z)dz = . (6) 

4 dz 4 

It is convenient to solve Eq. 10 for the local variation in volume fraction, Aip(z) — 
if(z) — ifo, rather than directly for tp(z). Assuming again linear elasticity and using the 
boundary condition Q, after some standard manipulations one finds the equilibrium 
concentration profile: 

z — h. 



Aip{z) = Vov 



1 — exp 



(7) 



where we have defined a third length £ related to the elastic and solid friction 
characteristic lengths introduced above by 1/1 = 1/L — 1/A. Note that there are 
no a priori restrictions on the sign of £; I = is also possible, if the elastic and 
friction characteristic lengths are equal. We will show in the following that for our 
gels < I < h n : in this case Eq. Q predicts the concentration profile to be essentially 
constant, except for a segment of length a few ts at the top of the gel. In general, 
we expect both A and L to depend on the gel volume fraction ipo; additionally, L 
may contain any deviations with respect to the usual Jannsen's law for ordinary solid 
materials, e.g. due to the fractal morphology of the gels. For L — > oo, solid friction 
becomes negligible and Eq. reduces to the exponential profile predicted by Eq. J5J) . 
By imposing particle conservation, A^ is found to satisfy 

A 00 = j{h 00 -e[l-eM-hoo/t)}}- (8) 

Because hoc — ho — Aoo, this expression is an implicit equation for Aoo when A and L 
are known and the initial height is assigned. By contrast, Eq. (JSJ) may be used directly 
to fit the experimental data, which can be easily plotted as Aoo vs hoc, rather than 
vs h as in Fig. 

Various settling regimes are possible, depending on the relative values of A, L, and 
ho (or, equivalently, hoc). Three limiting cases are identified, according to the behavior 
of the argument of the exponential term in the r.h.s of Eq. 0: i) hoo/i — > 0, for which 
Aoo w h 2 OQ /(2\); ii) h^/l -> oo, for which Aoo ~ h^l/X- £ 2 /X; iii) h^jl -> -oo, for 
which Aoo ~ (A 2 /^) exp(/ioo/|^|). The assumption of linear elasticity that was made 
in deriving Eqs. Q and (jHJ poses an additional constraint on the solutions for Aoo, 
since linear elasticity only applies to small deformations Aoo ^ ^oo- We recapitulate 
the different settling regimes in Fig. 0] which shows schematically a "state diagram" 
of the solutions to Eq. (JSJ, using as "state variables" X/hca and L/hoo and classifying 
the various zones according to the value of L relative to A. Zone I corresponds to 
comparable friction and elasticity characteristic lengths (|A — L\ << XL/h^). In 
zones Ha and lib the elastic length scale dominates over the friction length scale 
{L/hoo << 1 and L/hoo « A/ft-oo for Ha and 1 << L/hoo « A//ioo for lib). The 
reverse applies to zones Ilia, where A/ft-oo << 1 and X/hoo << L/hoo, and 111b, for 
which 1 << X/hoo << L/hoo- In summary, for the L/hoo > 1 portion of zone I and for 
zones lib and 111b one finds that a quadratic law applies: Aoo ~ ^oo/(2A) ~ h 2 /(2X). 
Note that this is the same expression as derived in the absence of solid friction. Hence, 
in these zones and to leading order solid friction has no influence on the asymptotic 
settling. At first sight, this result is somehow surprising for zone lib, where L < X and 
one could have guessed that solid friction would dominate over elasticity. However, 
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Figure 4. "State diagram" illustrating the location of the various asymptotic 
settling regimes discussed in the text in the (logA/hoo, log L/hco) plane. The 
corresponding approximate behavior of Aoo as a function of hoo is shown on 
the right. The A < h^o portion of zone I and zone Ilia should be disregarded, 
since they correspond to a large deformation regime that is incompatible with the 
assumption of linear elasticity. 



our analysis shows that solid friction is negligible as long as L is much larger than the 
asymptotic height of the gel, regardless of the characteristic length for elasticity. In 
zone Ha, a linear growth is found: 

I t 2 

Aoo « ^oo -j, (9) 

or, equivalently, 

L L 2 

Ao. « -ho - j-j- . (10) 

As one can easily check, the second term of the r.h.s. of Eq. (|10|l is negative and its 
magnitude is small compared to h^L/X, yielding a nearly linear dependence of A M on 
ho, with a slightly negative ho — intercept. This is precisely the behavior observed 
experimentally for our gels (see Fig.[3]and related discussion), strongly suggesting that 
our samples belong to zone Ha, where both elasticity and solid friction are present, 
but the contribution of the latter dominates. Finally, the L/h^ < 1 portion of zone I 
and zone Ilia correspond to large deformations Aoo > hoo that are incompatible with 
the assumption of linear elasticity and thus should be disregarded. 

To further test our model, we extract the two characteristic lengths L and A from 
the height dependence of Aoo and analyze both their (fo dependence and absolute 
magnitude [25| . Their volume fraction dependence is shown in Fig. together with 
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Figure 5. Volume fraction dependence of a) the elastic characteristic length A, 
b) the solid friction characteristic length L, c) the ratio L/X. The lines are power 
law fits to the data, whose exponent is shown in the panels. The values shown 
here confirm that our gels belong to zone Ha of Fig. [I] as suggested by the linear 
dependence of vs ho. 



their ratio. The elastic length A increases as </>q' 7±0 ' 4 , in excellent agreement with 
the scaling reported for the shear modulus of similar colloidal gels ^3 HOI 123] • 
Moreover, for ipo = 10~ 2 we find A = 2.35 cm, from which we calculate C = 2.8 Pa. 
This is of the same order of magnitude of, albeit somehow smaller than, the clastic 
shear modulus G ss 7 Pa that we estimate by extrapolating the data of ref. t° 
cpa = 10~ 2 . The solid friction length scale L is shown in Fig. Eb). If one assumes 
that the relation L = D/{Aji s K) derived for granular materials by Jennsen may 
apply to the gels and furthermore uses fi s « 0.6 and K < 1 as found typically for 
macroscopic materials, one finds L > 0.7 cm, not much larger than L m 0.1 cm, the 
experimental value for the largest tpo. Remarkably, however, we find that L decreases 
with decreasing volume fraction, strongly suggesting that the fractal nature of the gel 
has to be taken into account to correctly describe solid friction. Note that the variation 
of L with tpo is somehow counterintuitive, since larger values of L correspond to lesser 
solid friction: hence, the higher the volume fraction, the smaller the solid friction the 
gel experiences. We will come back to this surprising behavior in the following. The 
ratio L/X is found to decrease with tp Q , a consequence of the strong increase of the 
elastic length scale with volume fraction. As a consistency test of the theory developed 
above, we check that the magnitude of A and L are compatible with the assumption 
that our gels belong to zone Ha of the state diagram shown in Fig. |H (£ << A and 
L « /ioo), as suggested by the linear dependence of Aoo vs ho. Indeed, we find 
that for all samples L/X < 0.54 and L/hoa < 0.12. Moreover, the small deformation 
requirement is reasonably well fulfilled, since < 0.33/io for all samples, with the 
exception of the most diluted one, for which Aoo = 0.54/irj- 

In order to explain the volume fraction dependence of the solid friction 
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characteristic length, we propose a simple argument based on the assumption that the 
solid friction stress between the wall and the gel is proportional to N, the number of 
particles per unit area in contact with the wall. During settling the gel is compressed; 
as a result some particles are pushed against the wall. The lower ipo, the weaker the 
gel and the larger the number of particles pushed on the wall, thus leading to enhanced 
solid friction and smaller values of L. This argument may seem in contradiction with 
Amonton's law which sates that the frictional force F is proportional to the load P 
and is independent of the area A of the surfaces in contact j2Sj- However, as first 
observed by Coulomb, in the presence of adhesion forces F is a linear function of 
both load and contact area: F = fiP + cA. As explained by Ringlcin and Robbins 
27 , if the real contact area A is proportional to the load, as it is the case for elastic 
solids, Amonton's law follows automatically: F = fj, s P with fi s = fi + kc. So, our 
assumption that the the solid friction is proportional to the number of solid particles 
per unit area in contact with the wall will hold if the following conditions are fulfilled: 
(i) the friction force between the gel and the container is mainly dominated by the 
adhesion of solid particles to the wall (kc 3> fi)', (ii) adhesive forces are negligible in 
determining the flattening of contact between the gel and the wall (the compaction of 
the gel is dominated by the gravity field). Both conditions seem realistic. When the 
gel is formed, the aggregates can be described as spherical blobs of volume V ~ £ 3 , 
with a negligible number of solid bonds in direct contact with the wall (zero initial 
friction) . As the equilibrium compaction of the gel is reached, the relative variation of 
the blob volume, which we assume to be isotropic for simplicity, can be approximated 
by the first term on the r.h.s. of Eq. (TOJ|: AV/V ~ A^/ho ~ L/X ~ A£/£, where the 
last relation applies to small deformations. For a fractal blob, the number of particles 
scales with blob size as n ~ £ df ; therefore, the number of solid bonds per blob which 
will touch the wall due to a change in size A£ is An ~ £ d -f _1 A£ ~ £ d fL/X. The 
number of new contacts per unit area is N ~ An/£ 2 ~ £ df ~ 2 L/X; using the scaling 
L/X ~ (fiQ 1 ' 5 obtained experimentally and £ ~ <^q > f° r 4f = 1-9 one finally finds 

L ~ N^ 1 ~ </?o' 4± °' 4 ) (11) 

where the uncertainty in the exponent has been calculated assuming an uncertainty of 
0.05 in the fractal dimension. This expression captures correctly the decrease of the 
solid friction between the wall and the gel as its volume fraction increases, as observed 
experimentally, although the exponent in Eq. Ijlljl is lower than the measured one. In 
view of the large experimental uncertainties in the scaling of L (see Fig. [SJ and the 
high sensitivity to the exact value of df of the exponent, the agreement seems however 
reasonable. 

As a final remak, we note that the model proposed in this work may explain 
some features in previously published data on the settling of strong gels. In Fig. 1 of 
reference [TOJ, Senis and Allain plot f2 = 1 — Aqo/Zio as a function of volume fraction. 
A power law behavior is observed, except for ip$ > 0.01, for which £1 tends to saturate 
to one, since for high- volume- fraction, stiff gels A^ << Iiq. One may wonder whether 
our model may apply to these concentrated samples. Indeed, we find that our data 
plotted as f2 vs ipo fit on the high volume fraction portion of the curve shown by Senis 
and Allain (except for an irrelevant shift factor in volume fraction, due to the difference 
in the physico-chemistry of the systems). This strongly suggests that for cpo > 0.01 
the gels of reference ^U] retain their integrity and can be described by our model. 
Furthermore, one can see in Fig. 2 of reference [TO] that for (p > 0.004 and h > 40 
mm m constant similarly to our samples, implying that Aqo oc ho and confirming 
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that the model proposed here applies also to the most concentrated samples of Senis 
and Allain. Similarly, we wonder whether some of the data by Manley and coworkers 
|13| may fit in the framework of our model. We recall that their system is very close 
to ours, the main difference being the salt that is used to induce aggregation (MgCl 2 
at 20 mM rather than NaCl at 2 M) and the range of volume fractions investigated 
(0.005 < ipo < 0.08 in as opposed to 0.0025 < ip < 0.01). The fact that most 
of their data at large volume fraction can be reproduced neglecting solid friction is 
consistent with our finding that L decreases with increasing ipo, making solid friction 
less important for more concentrated gels. However, for the less concentrated gel and 
the tallest samples the ho dependence of r shown in Fig. 3 of reference J3| clearly 
departs from the quadratic law predicted in the absence of solid friction. If a similar 
deviation was observed also for Aoo vs Jiq, this would suggest that for their most 
diluted and tallest samples solid friction does indeed play a role. 

5. Conclusions 

We have investigated the settling of rather diluted colloidal gels made of strongly 
bound silica particles, focussing on the small deformation regime where the network 
integrity is preserved. Contrary to previous works, we find that the asymptotic settling 
is incompatible with a model based only on a linear elastic response to the gravitational 
stress. By contrast, we are able to model our results by adding a solid friction term to 
the balance of forces acting on the gel. We find that solid friction plays a crucial role 
as long as the associated screening length L is smaller than both the gel height and the 
characteristic length scale A for elastic deformation. Therefore, solid friction should 
matter for tall gels at moderate to small volume fractions. By contrast, solid friction 
is negligible when L is comparable to or larger than the gel height, regardless of the 
magnitude of A. Additionally, we have shown that, due to the fractal morphology 
of the aggregates composing the gel, the contribution of solid friction becomes more 
important as the gel volume fraction decreases, a somehow surprising result. 

Our work provides a more general framework for understanding and modeling the 
asymptotic behavior of settling gels and could easily be extended to other viscoelastic 
media. Future investigations should address the full kinetics of sedimentation, where 
we expect solid friction to play again a role, together with viscous dissipation and 
elastic forces. 
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